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The luminosity function of high-redshift QSOs 
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Abstract. We measure the luminosity function of QSOs in the redshift range 
3.5 < z < 5.2 for the absolute magnitude interval —21 < -M145 < —28. Suit- 
able criteria are defined to select faint QSOs in the GOODS fields, checking their 
effectiveness and completeness in detail. The confirmed sample of faint QSOs 
is compared with a brighter one derived from the SDSS. Using a Monte-Carlo 
technique we estimate the properties of the luminosity function. Our results 
show that models based on pure density evolution show better agreement with 
observation than models based on pure luminosity evolution, even if a different 
break magnitude with respect to z ~ 2.1 is required at 3.5 < z < 5.2. Ac- 
cording to our modeling a faint-end slope steeper than low-redshift observations 
is required to reproduce the data, moreover models with a steep bright-end slope 
score a higher probability than models with a bright-end flattening. Determin- 
ing the faint-end of the luminosity function at these redshifts provides important 
constraints on models of the joint evolution of galaxies and AGNs. 



Introduction 



In recent years evidence has grown that QSOs represent a necessary phase of 
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Monaco. Fontanot & Taffonil l2006h are able to repro- 



duce the properties of present-day elliptical galaxies by assuming that their vig- 
orous star formation at high redshift is quenched by the feedback (i.e. galactic 
winds). This process is likely triggered by the QSO activity fed by the accretion 
of matter ont o a super massive black ho le (SMBH) at the center of the galaxy 
itself (see also Monaco Sz Fontanot 20051 ). In this framework the detailed deter- 
mination of properties of high redshift QSOs is fundamental in understanding 
the phenomena driving galaxy formation. In particular the study of the Lumi- 
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nosity Function (LF) provides strong constraints for the physical mechanisms 
involved. 

A great improvement in the determination of this statistical quantity has 
been achieved thanks to the large amount of data coming from major obser- 
vational programs such as the Two Degree Field QSO Redshift Survey (2QZ, 
ICroom et al. 20041 ) and in particu lar the third edition o f the Sloan Digital Sky 
Survey Quasar Catalog (DR3QSO. ISchneider et aJ]|2005l ) . However, at high red- 
shift the SDSS is sensitive only to the brightest QSOs; on the other hand the 
faint-end of the high-z QSO LF plays a key role for comparing different predic- 
tions of the formation and evolution of galaxies. In this work we estimate the 
QSO LF and space density at 3.5 < z < 5.2, down to M145 < —21. 



The database 



We combine a bright quasar s ample extracted from the DR3QSO, with a faint 
sample recently determined bv lCristiani et all ( 20041 ) and lFontanot et all ( 2006al ) 
using optical and X-ray observations in the framework of the Great Observato- 
ries Origins Deep Survey (GOODS) project. The selection of QSO candidates 
at different redshift in SDSS field has been carried out using color criteria: 
Richards et al.l (2002) presented the updated version of these criteria and in 
the following we use their eq. 6 and 7 as SDSS selection criteria to select a 
high-z QSOs sample out of the DR3QSO. The DR3QSO is a sample with an 
incomplete spectroscopic follow-up. We consider those objects in the SDSS pho- 
tometric catalogue that satisfy the SDSS criteria and we compare them to the 
corresponding SDSS spectroscopic catalogue: this way, we estimate that ~ 21% 
of the high-z QSO candidates have a spectrum. 

In order to select a sample of faint QSOs from the GOODS catalogues 
we consider the optical catalogues obtained usin g the ACS (Advanced C amera 



for Surveys) in the -B445, V66o> 2775, z 850 bands ( Giavalisco et alj|2004al ). The 



catalogues are z-band-based with magnitude limits in the four bands are re- 
spectively 27.50, 27.25, 27.00, 26.5. The selection of the QSO candidates is 
carried out in the magnitude interval 22.25 < zgso < 25.25. Expected QSO 
colors in different photometri c systems are es t imate d as a function of redshift 
using a template spectrum of ICristiani Viol ( 19901 ) for the QSO spectral en- 
ergy distribution (SED) convolved with a model of the intergalactic medium 
(IGM) absorption. Four optical criteria has been tailo red to select QSOs in 



the redshift interva l 3.5 < z <^ 5.2 as described in ICristiani et al.l (|2004l ) 
and iFontanot et al. ( 2006a ). These color selection criteria are known to se- 
lect a broad range of high-z AGN, not limited to broad-lined (type-1) QSOs, 
and are less stringent t han those typically used to identify high-z galaxies (e.g. 
Giavalisco et al. 2004bl ) . In order to reduce the number of The optically selected 
candidates we match them with X-ray C handra observations in t he HDF-N and 
CDF-S (for 2 Ms and 1 Ms, respectively lAlexander et al.l [20031 : iGiacconi et al. 
2002), within an error radius corresponding to the 3 <x X-ray positional uncer- 
tainty. The X-ray completeness limits over ~ 90% of the area of the GOODS 
fields are similar, with flux limits (S/N= 5) of ~ 1.7x 10 -16 erg cm -2 s — 1 (0.5- 

" erg cm -2 s _1 (2-8 keV) in the HDF-N field and 
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(2-8 keV) in the CDF-S field. The sensitivity at the aim point is about 2 and 
4 times better for the CDF-S and HDF-N, respectively. Given those flux lim- 
its, Type-1 QSOs with M145 < —21 are detecta ble up to z > 5.2, up to an 
optical-to-X-ray flux ratio q o;c ^ —1.7. Recently ISteffen et al.l ( 20061 ) provide 
an estimate of a ox statistics (their Table 5): at this optical luminosity they 
observed a mean value lower than this limit (a ox = —1.408 ± 0.165). Con- 
versely, any z > 3.5 source in the GOODS region detected in the X-rays must 
harbor an AGN (-M0.5 - 2 keV) ^ 10 43 erg s" 1 ). The final sample of QSO 
candidates c onsists of 16 candidate s, 10 in CDF-S and 6 in HDF-N and it is 
presented in lFontanot et al.l ( 2006al . their Table 1). 



S pectroscopic information for all 16 cand i dates has been obta i ned jV anzclla et al 
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Cristiani et al.l l2000l ) . Thirteen (5 in the HDF-N and 8 in CDF-S) can 
didates turn out to be AGN, with 8 QSOs at z > 3.5 (3 and 5), of which 
2 (0 and 2) are identified as Type II. One object turns out to be a QSO with 
a redshift z == 4 .759. In addition to the already known QSO at z = 5.189 
( Barger et al.l I2OO2I ) . it brings to two the total number of objects in the red- 
shift range 4 < z < 5.2, where our selection criteria are expected to be most 
complete and reliable. 



3. Estimating the luminosity function of high-z QSOs 

In order to build up the high-z QSO LF from the joint analysis of GOODS and 
SDSS observations a Monte-Carlo technique is adopted. Mock QSO catalogues 
in GOODS and SDSS photometric systems are generated, starting from an as- 
sumed parameterization of the LF and they are compared to real catalogues 
using a x 2 estimator to quantify the agreement between real data and the func- 
tional form of the LF. In order to clone the properties of high-z objects we define 
a QSO template library, tailored using high-quality SDSS QSO spectra at lower 
redshift. We choose the redshift interval 2.2 < z < 2.25, for which the SDSS 
sample had the highest possible level of completeness and the continuum of the 
QSOs was sampled in the largest possible wavelength interval from the Ly a 
emission upward. The rest-frame spectrum of each selected object blueward of 
the Ly n lin e is es timated using the continuum fitting technique derived from 
Natali et al.l ( 19981 ). The library is then used to compute the distribution of the- 



oretical QSOs colors at increasing redshift in the SDSS and ACS photometric 
systems. In Fig. Q] (left panel) we show our results for the SDSS photometric 
system compared to the colors of observed QSOs in DR3QSO. The points refer 
to the colors of observed QSOs at different redshifts (red stars refer to objects 
with 3.5 < z < 4.0; green squares to 4.0 < z < 5.2; filled circles to z > 5.2; 
cyan dots to z < 3.5). The dotted line shows the projection of the SDSS selec- 
tion criteria. The mean QSO color in our library is represented by the solid line. 
Dashed lines are 5% and 95% percentiles of color distribution in the template 
library (photometric errors are not accounted for). Similar results hold for the 
GOODS survey. 



In order to estimate the QSO LF, the approach of lLa Franca &: Cristiani 



( 19971 ) is adopted. An LF of the form of a double power law is assumed, a pure 



luminosity evolution (PLE) or a pure density evolution (PDE) with a power 
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Figure 1. Left panel: Color diagrams for the confirmed QSOs in DR3QSO. 
Right panel: Completeness of selection criteria at various redshifts. The solid 
line refe rs to the SPSS selecti on criteria for SDSS sources. The dashed line 
refers to ICristiani et al.l (|2004f ) selection criteria in the GOODS fields. 



law form and, alternatively, with an exponential form are adopted. Given a 
value for the break luminosity M*, the slopes a and f3 of the double power law, 
the normalization and the redshift evolution parameter k z , it is possible to 
calculate the expected number of objects from the LF up to a given magnitude 
in a given area of the sky. For each object a value of absolute M145, a redshift 
and a template spectrum are extracted. We then use the template k-correction 
and colors at the corresponding redshift to simulate the apparent magnitudes 
in the SDSS and ACS photometric systems. The photometric errors in each 
band are estimated using SDSS and GOODS main photometric catalogues, and 
added to the simulated magnitudes. The corresponding color selection criteria 
are applied to the generated catalogues to obtain mock SDSS and GOODS 
QSO catalogues. In order to limit the statistical error multiple realizations 
are considered for each parameter set. The agreement between real and mock 
catalogues is evaluated using both a % 2 statistic and a Kolmogorov-Smirnoff 
test. The best fit parameters for each model are estimated usi ng a minimization 
technique. The results are collected in lFontanot et all ( 2006a! ) . Table 3 and the 
best-fitting model is shown in Fig. [2] (left panel) . 



4. Results 



It is impossible to reproduce the observed QSO distribution with PLE models. 
PDE models show better agreement with observations, even if a luminosity 
evolution of the magnitude of the break, M^ 45 , from z = 2.1 to z — 3.5 is 
required. Parameterizati ons assuming a relati vely steep faint-end slope similar 
to the low-z estimate by iRichards et all ( 20051 . SDSS combined to 2dF sample) 
score a higher probability with respect to those with a faint-end slope as flat 



as m 



Croom et al. ( 2004L 2dF sampl e). Models wit h a shallower slope of the 



bright-end of the LF (as suggested bv lFan et al.ll2003l ) are not able to reproduce 
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Figure 2. Left panel: analytical fit to the high-z QSO LF. The solid line 
shows, as a reference, the best fit model. The position of the filled sym- 
bols is obtained by multiplying the values of this model by the ratio be- 
tween the number of observe d sources and the number of simulated sources. 
lLa Franca fc Crist iani (1997|) demonstrated that this technique is less prone 
to evolutionary biases with respect to the conventional 1/V max technique. 
Right panel: Spac e density evolution of QSOs, compared to the predictions o f 
MORGANA model (|Monaco. Fontanot & Taffonill2006t iFontanot et al"1l2006bf) . 
The (black) solid line, the (magenta) dotted line, the (blue) dashed line and 
the (red) dot-dashed refer to different values of the <r parameter (see text): 
0, 30, 60 and 90 respectively (values are in Kms -1 ). 



the observed distributions: a bright-end sl ope as steep as i n Croom et al.l ( 20041 ) 
is required. This result is at variance with lRichards et al.l (|2006l h 

In order to understand the origin of this discrepancy we focus our attention 
on the completeness of the samples selected using the GOODS and SDSS criteria. 
To obtain a robust estimate of this quantity we apply the selection criteria to our 
template library and analyze the fraction of QSOs recovered at various redshifts. 
The resulting c o mplet eness is shown in Figure [T] (right panel) . As a comparison 
iRichards et al.l (200(3) estimated the completeness of their criteria to be well 
above 90% in the whole range of redshift of interest. We check that this is related 
to the fact that the SDSS selection criteria are tailored on a QSO template 
spectra whose mean continuum slope f u oc v f is "bluer" (7 = 0.5 ± 0.3 



Fan et al.|[l999l ) than the mean slope in our template library (7 = 0.7 ± 0.3 



The inferred completeness has a direct co nsequence on th e shap e and evolution 
of the estimated QSO LF. Assuming the IRichards et al.l ( 20061 ) completeness, 
the models with a shallow bright-end score higher probabilities with respect to 
models with a steep bright-end. 

We use our best estimate for the LF to compute the QSO contribution to 
th e UV backgro u nd at redshift 3.5 < z < 5.2, following the method outlined 
bv lBarger et al. ( 20031 ). We conclude that the QSO contribution is insufficient 
for ionizing the IGM at these redshifts. 

When compared to physically motivated models, the present results show 
that to reproduce the observations it is necessary to suppress the formation/feeding 
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of low- mass BHs inside DMHs at high redshift. As a consequence the present es- 
timate of the space density of high-z faint QSOs gives strong constraints to mod- 
els of joint formation and evolution of galaxies and AGN. In fig. Q] (right panel) we 
compare the evolution of the QSO space den sity in bin of bolometric luminosity 
with the predictions of the MORGANA model ( Monaco. Fontanot Taffonill2006l . 
the detailed description of the model i s beyond the scope of pr esent proceedings 
and can be found in that paper). In iFontanot et all ( 2006bl ) we propose that 
kinetic feedback is the key process responsible for the downsizing of the AGN 
population: in that paper we demonstrate that the velocity dispersion of cold 
clouds <r, due to turbulence induced by exploding SNe, scales with the star for- 
mation time-scale i* as cr = <xo(i*/l Gyr) 1 / 3 ; the normalization ctq regulates 
the level of turbulence, so that a high value (~ 50 km s _1 ) induces significant 
mass loss from small star-forming bulges at high redshift. Different lines in fig.Q] 
(right panel) refer to model with different values for <xo (see caption for more 
details). It is evident from the figure that the space density evolution of high- 
z low-luminosity AGN is of fundamental importance to constrain this physical 
mechanism. 
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